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The macroscopic mechanical properties of densely packed coagulated colloidal particle gels 
strongly depend on the local arrangement of the powder particles on length scales of a few par- 
ticle diameters. Heterogeneous microstructures exhibit up to one order of magnitude higher elastic 
properties and yield strengths than their homogeneous counterparts. The microstructures of these 
gels are analyzed by the straight path method quantifying quasi-linear particle arrangements of par- 
ticles. They show similar characteristics than force chains bearing the mechanical load in granular 
material. Applying this concept to gels revealed that heterogeneous colloidal microstructures show a 
significantly higher straight paths density and exhibit longer straight paths than their homogeneous 
counterparts. 
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I. INTRODUCTION 

Mechanical properties of coagulated colloids are im- 
portant in many technical areas. Sediments [lj, ceramic 
pastes and suspensions [2], pharmaceutical formulations 
such as cremes and emulsions Q and some food [|| are 
examples, for which it is desirable to control the mechan- 
ical properties. 

Recently, an internal gelation method (DCC = di- 
rect coagulation casting) [j| Q was developed to pro- 
cess electrostatically stabilized colloidal suspensions to 
coagulated particle gels. Thereby it was found that this 
method permits to control the gels' microstructures for 
volume fractions ranging from 0.2 to 0.6. The method 
allows for an in situ, i.e., undisturbed destabilization of 
the colloidal suspension by either changing the pH of the 





FIG. 1: Microstructure of coagulated silica suspensions at 
a volume fraction of 0.4 formed by the ApH- (left) and the 
Al-method (right) using the DCC process (particle diame- 
ter 0.525 fim). Micrographs obtained by cryogenic scanning 
electron microscopy [7(. 
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solution (ApH-method) resulting in a "homogeneous" 
microstructure or by increasing the ionic strength (AI- 
method) leading to a "heterogeneous" microstructure 0] • 

Figure [T] shows cryogenic scanning electron microscopy 
(cryo-SEM) pictures of both microstructures. The parti- 
cles in the ApH-coagulated gel present a highly homoge- 
neous microstructure whereas in the AFsystem inhomo- 
geneities on the length scale of a few particle diameters 
are observed. This qualitative observation was quantified 
by the three-dimensional pair-correlation function calcu- 
lated from stereo cryo-SEM images 0] and by in situ per- 
formed diffusing wave spectroscopy (DWS) experiments 
during destabilization [8| . The structural differences cor- 
respond to differences in the heterogeneity: the first peak 
of the pair-correlation function at r/d = 1 (r being the 
distance between the particles and d the particle diame- 
ter) is higher for the Al-system indicating locally denser 
regions with higher average coordination number. Also, 
local maxima are present in the pair-correlation function 
for the Al-system at r/d » 1.4 and at r/d » 1.6 corre- 
sponding to characteristic peaks in hexagonally packed 
particle arrangements. During the Al-destabilization, 
particle rearrangements lead to larger pores and thus 
to less homogeneous microstructures [8|, whereas during 
the ApH-destabilization, the initially stabilized, liquid- 
like microstructure, is "frozen" , which results in a more 
homogeneous microstructure. 

Rheological and uniaxial compression experiments on 
coagulated colloidal structures obtained by DCC @, [Toj 
showed that colloids with heterogeneous microstruc- 
tures have significantly higher elastic moduli and yield 
strengths than their homogeneous counterparts. Mea- 
suring the dynamics of the colloids during destabilization 
by DWS confirmed these findings quantitatively using a 
model for the storage modulus proposed by Krall and 

Weitz d, [Trj]. 

Alkali-swellable polymer (ASP) particles were used 
in order to introduce heterogeneities during the pH- 
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FIG. 2: Elastic plateau storage modulus G' p of alumina par- 
ticle suspensions (average particle diameter do = 0.4 /jm) in 
dependence of the volume fraction formed by the ApH- and 
the Al-method of the DCC process [Iol |. 

destabilization. Small amounts of 80 nm ASP particles 
were admixed to the powders under acidic conditions. 
These particles swell upon changing pH during the in- 
ternal gelling reaction of the DCC process and enfold 
to 0.7 /xm size producing less homogeneous microstruc- 
tures. Those samples with swollen polymer particles 
showed much higher mechanical properties in compar- 
ison to samples without swellable polymers and hence 
more homogenous microstructures [12j . In particular, 
they present the same high mechanical properties as sam- 
ples with heterogeneous microstructures produced by the 
Al-method [H[. 

Examples of the rheologically measured elastic proper- 
ties of alumina particle suspensions (average particle di- 
ameter do = 0.4 (an) for different volume fractions desta- 
bilized by the ApH- and the Al-method, respectively, are 
shown in Fig. [2] Almost one order of magnitude higher 
elastic plateau storage moduli are measured for hetero- 
geneous microstructures than for those with homogenous 
microstructures [To| . 

In summary, strong evidence is given that the dif- 
ferences in macroscopic mechanical properties of coag- 
ulated particle suspensions are controlled by the differ- 
ences in microstructure. An open question is now how 
these microstructural differences on particle length scales 
can have such a dramatic influence on the mechanical 
properties. Therefore, a concise quantitative analysis of 
microstructures of colloidal particle systems is needed. 

In preceding works, various characterization methods, 
such as the radial pair-correlation function [l3j , the bond 
angle distribution function [l3| . the triangle distribution 
function [13l ] and the Minkowski functionals in conjunc- 
tion with the parallel-body technique [l4T| . were applied 
to sets of microstructures generated by Brownian dynam- 
ics (BD) simulations [l5j]. These simulations were used 
to study the coagulation dynamics and the evolving mi- 
crostructures in dense colloidal suspensions and the re- 



sulting microstructures agree well with experiments [16| . 

While the pair-correlation function permits to quan- 
tify the amount of structural rearrangement during the 
coagulation the analysis using the Minkowski func- 
tionals in conjunction with the parallel-body technique 
supplies additional information on the structure's mor- 
phology resolving microstructural differences on a length 
scale limited by the largest pore size flij . The bond 
angle distribution function and the triangle distribution 
functions are useful means to examine the local building 
blocks of the particle network [l3j]. Particular features, 
as, for example, peaks in the respective distribution func- 
tion have successfully been correlated to the structure's 
heterogeneity and porosity. The same conclusion is valid 
for the Minkowski functionals [Til ]. 

All these four methods are good means to compare 
structures in terms of their heterogeneity. However, 
these structural descriptions do not unambiguously help 
to understand why more heterogeneous colloidal struc- 
tures possess stronger mechanical properties, as they do 
not adequately capture the microstructural characteris- 
tics that are responsible for the mechanical properties of 
these particle systems. 

It is well known from granular matter physics that par- 
ticulate systems under mechanical stress carry load via 
chains of contacting particles, termed as force chains, as 
observed in ex per iments on granular materials [l7| and 
in simulations [18[. Granular materials, as, for example, 
sand piles, that contain large amounts of particles ar- 
ranged in chains of contacting particles are expected to 
possess higher mechanical properties than those in which 
the particles have to rearrange upon applied external load 
in order to form such chains. Hence it is plausible that 
these chains of contacting particles control the mechani- 
cal properties. 

The aim of this paper is to adequately analyze the 
microstructure of dense colloids and to link microstruc- 
tural differences to differences in macroscopic mechanical 
properties upon applied mechanical stress. Thereby, the 
focus is on the identification of densely packed regions 
and chains of contacting particles in these particle net- 
works that may correlate differences in mechanical prop- 
erties with differences in microstructures. Our approach 
is to analyze the distribution of densely packed regions in 
each microstructure using the common neighbor distribu- 
tion and the dihedral angle distribution function. Addi- 
tionally, a new method called the straight path method, 
characterizing the quasi-linear arrangement of particle 
chains, is introduced. The evaluation of these character- 
ization methods is performed with regard to their abil- 
ity to distinguish quantitatively between homogeneous 
and heterogeneous microstructures obtained from BD- 
simulations [l5| mimicking well those found in experi- 
ments [l6| . 
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TABLE I: Potential Parameters for the Coagulating Suspen- 
sions 



Parameter 




Symbol Value 


Hamaker constant of AI2O3 


in water 


Ah 


4.76 x 10" 20 J 


Particle diameter 




d 


0.5 /jm 


Relative dielectric constant 


of water 


tr 


81 


Absolute temperature 




T 


293 K 


Valency of ions 




Zv 


1 


Inverse Debye screening len 


5th 


n 


10 8 m _1 



II. MATERIALS AND METHODS 
A. Structure Generation 

The homogeneous and the heterogeneous microstruc- 
tures are fully coagulated colloidal suspensions obtained 
from Brownian dynamics (BD) simulations [l5j . The 
Derjaguin-Landau-Verweg-Overbeek (DLVO) theory (I9I ] 
was used to describe the particle-particle interaction, 
where pair-wise particle potential interactions are as- 
sumed given by the sum of the van der Waals attraction 
yvdw ^gq^ jjjj) anc j t ne electrostatic double layer repul- 
sion V el (Eq. ©). Thus, V dlvo = V vdw + V el with 
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FIG. 3: DLVO interaction for different surface potentials: 
(a) pure van der Waals attraction (ApH-destabilization), (b) 
AE = 0fc fl T and (c) AE = 5.65fc s T (Al-destabilization). 
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(1) FIG. 4: Configuration (grey particles) with two common 
neighbors (white particles) separated by a distance r and hav- 
ing a dihedral angle £. Perspective view (a) and top view (b). 



V el (r) =ne r e 



tanh 



4hT 



(2) 

respectively. The DLVO parameters are summarized in 
Table |T] and the potential curves for various values of the 
surface potential ^0 are shown in Fig. [3J 

For = mV, the electrostatic double layer repulsion 
is zero and the inter-particle potential is only given by the 
attractive van der Waals potential. For ^0 = 12 mV a 
secondary minimum appears and AE denotes the energy 
barrier between the local maximum and the secondary 
minimum. For ^0 — 15 mV, a repulsive barrier of AE 
= 5.65 fcgT exists and the secondary minimum is found 
at r = 1.08d. The model further contains the frictional 
Stokes' drag force and a random Brownian force caused 
by the suspending liquid. 

We analyze microstructures generated using = 
mV, which correspond to the ApH destabilization 
method and ^0 = 15 mV, corresponding to the AI 
destabilization method. In the following, the ApH- 
microstructures are labeled as homogeneous and the AI- 
microstructures as heterogeneous. All structures have 



! the same volume fraction of 0.4 and contain 8000 spher- 

dexp (— n{r — d}) , ical particles forming one percolating cluster. The parti- 
cles have a diameter of 0.5 /xm. 



B. Structure Characterization Methods 



1. Common Neighbor and Dihedral Angle Distribution 
Function 



The common neighbor analysis [20j considers pairs of 
particles, referred to as configurations, and determines 
the number of particles that are in contact with both par- 
ticles of the configuration. In two dimensions, configura- 
tions can have at most two common neighbors. In three 
dimensions, configurations can have up to five common 
neighbors. The number of configurations with n com- 
mon neighbors is denoted CN n with n — 1, ... ,5. In this 
study, absolute numbers of configurations are compared 
as the microstructures analyzed have identical volume 
fraction, the same number of particles and equal particle 
diameters. 
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The common neighbor distribution describes the short 
range arrangement of particles. In particular, CN is the 
number of pairs of contacting particles with no common 
neighbor. CN\ counts the number of configurations with 
exactly one neighbor common to both particles, forming 
equilateral triangles. For n — 2 the four particles of the 
configuration form a regular tetrahedron if the two com- 
mon neighbors are in contact, and a generalized tetrahe- 
dron, having one longer edge, otherwise. 

Configurations with two or more common neighbors 
(CN2, . . . , CN5) can be further characterized by the di- 
hedral angle distribution. This distribution analyses the 
arrangement of the neighboring particles around the two 
particles of the configuration, which form regular trian- 
gles with each of their common neighbors. For n > 2 
there are at least two common neighbors. In this case, 
the dihedral angle £ is defined as the angle between two 
planes spanned by triangles belonging to the same config- 
uration as is shown in Fig. [5] £ is measured between one 
triangle and the triangles formed with each other com- 
mon neighbor of the original couple of particles. Chem- 
istry uses this method to characterize molecular struc- 
tures 21]. In the field of granular matter, the dihedral 
angle distribution has been used to study the structural 
organization and correlations in very large packings of 
monodispersed spherical particles [22|, [23| . 

For two common neighbors in contact with each other, 
the four particles form a regular tetrahedron and the di- 
hedral angle is 70.5°. This angle constitutes the lower 
limit of the dihedral angle distribution. The upper limit 
is 180° when the four particles are in a plane and form a 
square. 

Formally, the dihedral angle £ between two common 
neighbors of a configuration having a separation distance 

r between each other is defined by £ = 2arcsin {^^J- 

The dihedral angle distribution is given by p(£ + <5£/2) = 
SN^, £ + <5£) where 6N^, £ + <5£) is the number of gen- 
eralized (open) tetrahedra with a dihedral angle between 
£ and £ + (5£ 



2. Straight Path Method 

A straight path is a quasi-linear chain of contacting 
particles. The shortest chain consists of three particles. 
In order for three particles to form a straight path they 
have to fulfill the following condition: The angle 8 be- 
tween the vector connecting the first two particles (Pi, 
P2) and the vector connecting the second two particles 
(P2, P3) must be smaller than a chosen threshold angle 
8 C . This is exemplified in Fig. [5j 

Analogously, a straight path of length I is one of length 
I — 1 that has another particle attached to one end and 
whose connection vector fulfills the condition to stay in 
a direction within a cone of angle 9 C from the vector 
connecting the previous two particles. The distribution 
of straight path lengths is obtained by determining the 




FIG. 5: Chain of three particles forming a straight path. 



number of straight paths SPi of length I = 3, . . . , l m ax, 
with l max the longest path under consideration. Absolute 
numbers are calculated and compared in the distribution 
of SPi. 

The angle 8 C determines the quasi-linearity and speci- 
fies the maximally allowed deflection from an absolute 
linear arrangement of the particles. The dependence 
of the straight path statistic on the angle 8 C is pre- 
sented as part of the results section and is used in or- 
der to choose a suitable angle 8 C . The distribution of 
straight paths for the homogeneous and heterogeneous 
microstructures is compared. The average path length 

SPjnean — 

T,i>3lS Pi /S/> 3 S Pi and the number of paths 
longer than a chosen length l Q is determined. The lat- 
ter is defined by SPi>i = S/>; S'P/ (lo = 4 and 5) and 
permits to quantify differences in the number of straight 
paths with longer lengths, neglecting paths shorter than 
lo particles. 



III. RESULTS AND DISCUSSION 

A. Common Neighbor Analysis 

The number of particle configurations CN n is higher 
for the heterogeneous than for the homogeneous mi- 
crostructure for n = 0, ... ,3 (Fig. [5]). This can easily 
be understood as the mean coordination number of a 
particle in the heterogeneous microstructure is with 5.2 
roughly 10% larger than in the homogeneous microstruc- 
tures with 4.7. When normalized by the mean coor- 
dination number, the distributions collapse, indicating 
that the relative distribution of common neighbor con- 
figurations is identical for both microstructures. Thus, 
the common neighbor distribution does not capture in 
a suitable way the structural differences between these 
microstructures. 
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FIG. 6: Common neighbor distribution for the homogeneous 
and heterogeneous microstructures: number of configurations 
CN n in dependence of the number of common neighbors n. 
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FIG. 7: Dihedral angle distribution for the homogeneous and 
heterogeneous microstructures. For clarity purposes, only ev- 
ery fourth data point is indicated by a symbol. 



B. Dihedral Angle Distribution 

The dihedral angle distributions of both microstruc- 
tures show a characteristic peak at 70.5° (Fig. [7]) cor- 
responding to configurations that form regular tetrahe- 
dra with two common neighbors. The peak height is 
the same for both microstructures. The heterogeneous 
microstructure has a peak at 77°, which is missing for 
the homogeneous one. This angle corresponds to a sepa- 
ration distance between the common neighbors of 1.08<i 
which is the distance between particles trapped in the 
secondary minimum (Fig.[3|). 

At £ = 109.5° the homogeneous microstructure has 



a peak that corresponds to a common neighbor separa- 
tion of 1.41<i ~ V2d, which is absent in the case of the 
heterogeneous microstructure. This separation distance 
corresponds to two configurations forming a square and 
sharing two common neighbors. An additional character- 
istic peak is found at 141° having roughly the same height 
for both microstructures and corresponding to two juxta- 
posed tetrahedra. The peak at 148° is solely present for 
the heterogeneous microstructure and corresponds to one 
regular tetrahedron juxtaposed to one with a distance of 
1.08d between the common neighbors. 

In conclusion, even if there are sizable differences be- 
tween homogeneous and heterogeneous structures, the 
absolute number of regular tetrahedra is virtually equal 
for both microstructures and the analysis does not reveal 
any structural difference that could explain the better 
mechanical properties of the heterogeneous microstruc- 
tures. Thus, the common neighbor and the dihedral an- 
gle distributions are not suited to link the microstructural 
differences to the different mechanical properties. 



C. Straight Path Analysis 

The straight path analysis depends on the threshold 
angle 9 Cl hence, a suitable value 9 C has to be determined. 
For small angles 9 C only a small number of straight paths 
are expected because only a few particle chains satisfy 
the criterion to be a straight path. For large angles 9 C 
many straight paths are found but the physical mean- 
ing of a straight path is lost as the path becomes less 
"straight" and thus loses its load bearing capacity. Also 
for increasing 9 C the ratio of the number of paths with 
equal lengths between the heterogeneous and the homo- 
geneous microstructure decreases. Figure [5] shows the 
number of straight path longer or equal to five parti- 
cles for the homogeneous (SP^f) and for the hetero- 
geneous microstructure (SP^f) as well as their ratio 



(SP^,f : SP^) m dependence of the cone angle 9 C 
ranging from 10° to 40°. The lines serve as guide for 
the eyes. 

For 9 C between 14° and 18° the ratio SPgf : SPg° 
shows a maximum. However, below 22° the number of 
straight paths is quite small (SPi>5 < 100) for both mi- 
crostructures. Toward larger angles the absolute number 
of straight paths is increasing while the ratio SP^f : 
^Pi>?5 i s decreasing. The angle 9 C = 30° was chosen 
as this angle gives reasonable values for both the ratio 
( SP i>5 '■ SP i>5 ~ 3 ) and absolute number of paths 
(5P|f = 451andSP|° = 163). 

For both microstructures the distribution of straight 
paths shows an exponential decrease with straight path 
length Is p (Fig. [5]) . Both microstructures have approxi- 
mately the same number of straight paths of length three. 
Toward longer path lengths the microstructures present 
a diverging behavior: the heterogeneous microstructure 
has more paths of longer length than the homogeneous 
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FIG. 8: Number of straight paths for both microstructures 
longer or equal to five (SP^f and SPfl^ , open symbols) and 
their ratio (SPiLs '■ SPiL® , triangles) as a function of the 
threshold angle 9 C . 
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FIG. 9: Number of straight paths SPi in dependence of the 
straight path length Isp for the homogeneous and heteroge- 
neous microstructures (squares and circles, respectively) and 
exponential fits (dashed and solid line, respectively). 



TABLE II: Average path length SP m ean and numbers of paths 
longer or equal to four (SPi>4) and five (SPi>$) for both 
microstructures 



Homogeneous Heterogeneous 



S Prnean 


3.28 


3.45 


SP>4 


779 


1385 


sp> 5 


163 


451 



heterogeneous microstructure longer chains SPi>4 and 
SPi>5 are found 1.8 and 2.8 times more often, respec- 
tively, than in the homogeneous one. 

To conclude, the number of straight paths in the homo- 
geneous and heterogeneous microstructures follows an ex- 
ponential distribution and a higher characteristic inverse 
path length A is found for the homogeneous microstruc- 
ture. The heterogeneous microstructure contains signif- 
icantly more straight paths of longer lengths than the 
homogeneous one: in absolute numbers twice as many 
paths with a length > 4 and three times as many having 
a length > 5 are observed. These are very significant 
differences between the two microstructures. 

Our results suggest that the differences in the straight 
path distributions are characteristic for the better me- 
chanical properties of the heterogeneous microstructures. 
Indeed, in granular materials load is transmitted via force 
chains that are quasi-linear substructures of the parti- 
cle network. This has been experimentally verified us- 
ing photo-elasticity techniques [l7| and computationally 
modelled using the discrete element method [l8]. Also, 
the computational study of the force chain network in 
2D granular assemblies of polydispersed particles sub- 
jected to indentation by a rigid flat punch has shown 
that the force chain length follows an exponential distri- 
bution [H. 



IV. CONCLUSIONS 



microstructure. The distributions were fitted by an ex- 
ponential law using SPi cx exp(— XI), with A the char- 
acteristic inverse path length. A remarkably good fit of 
the data is obtained for both microstructures indicated 
by the values close to one for the correlation coefficients 
R 2 . The higher exponent A found for the homogeneous 
microstructure demonstrates the more rapid decrease of 
the number of straight paths toward longer path lengths 
than in the case of the heterogeneous microstructure. 

The values for the average path lengths (SP mean ) and 
the numbers of paths longer or equal to four (SPi>4,) 
and five (S , P/>s) are summarized in Table |ll] The av- 
erage path length for the heterogeneous microstructure 
is only 5% longer than for the homogeneous microstruc- 
ture because SPi=3 dominates the statistics and is almost 
identical for both structures. On the other hand, in the 



Different chemical pathways during the internal desta- 
bilization of colloidal suspensions lead to wet ceramic 
green bodies with drastically different mechanical prop- 
erties. The microstructures of these destabilized colloids 
reveal differences in the heterogeneity of the particle ar- 
rangements on a length scale in the order of a few parti- 
cle diameters. Various microstructural characterization 
methods, i.e., the common neighbor distribution, the di- 
hedral angle distribution and the straight path method, 
have been applied in order to analyze homogeneous and 
heterogeneous microstructures generated by Brownian 
dynamics simulations, which were shown to agree well 
with the experimentally determined microstructures of 
such coagulated colloidal particle systems. 

Toward our goal, to establish a correlation between the 
microstructural differences and the differences in macro- 
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scopic mechanical properties, it was found that the com- 
mon neighbor and the dihedral angle distributions do 
not discriminate enough the differences between the ho- 
mogeneous and heterogeneous microstructures and can 
therefore hardly account for the large differences in me- 
chanical properties. The newly introduced straight path 
method reveals significantly more straight paths of longer 
lengths for the heterogeneous than for the homogenous 
microstructure: twice as many straight paths of length 
> 4 particles and three times as many of length > 5 
particles are found in the heterogeneous microstructure. 
These differences in straight path number and length 
seem suitable to characterize and to differentiate between 
these structures and these differences are considered to be 
sufficiently large to account for the differences in mechan- 
ical properties. The straight paths may capture best the 
characteristic microstructural features which are relevant 
for the mechanical properties. The number of straight 
paths was found to follow an exponential distribution just 
like the distribution of the force chain lengths in mechan- 
ically loaded granular material [24]. Additionally, the 



quasi-linear structure of the straight paths seems to cor- 
respond to the geometrical shape of force chains [13, EH 
that are well known to determine the load bearing capac- 
ity of granular matter. These findings are encouraging for 
an attempt to establish a correlation between defined mi- 
crostructural features and the differences in mechanical 
properties of destabilized colloids. 

Further quantitative structural analyses and computa- 
tional modelling using the discrete element method are 
in progress to study the force chains in ceramic green 
bodies and to correlate them with the straight paths of 
their microstructures. 
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